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Abstract 

The s/g(2)-quantum group invariant spin 1/2 XXZ-Heisenberg model with open 
boundary conditions is investigated by means of the Bethe ansatz. As is weh known, 
quantum groups for q equal to a root of unity possess a finite number of "good" rep- 
resentations with non-zero q-dimension and "bad" ones with vanishing q-dimension. 
Correspondingly, the state space of an invariant Heisenberg chain decomposes into 
"good" and "bad" states. A "good" state may be described by a path of only "good" 
representations. It is shown that the "good" states are given by all "good" Bethe 
ansatz solutions with roots restricted to the first periodicity strip, i.e. only positive 
parity strings (in the language of Takahashi) are allowed. Applying Bethe's string 
counting technique completeness of the "good" Bethe states is proven, i.e. the same 
number of states is found as the number of all restricted path's on the sZg(2)-Bratteli 
diagram. It is the first time that a "completeness" proof for an anisotropic quantum 
invariant reduced Heisenberg model is performed. 



1 Introduction 

The Bethe ansatz method has been applied to a large number of integrable models as one- 
dimensional Heisenberg spin chains and statistical lattice models in two dimensions. The 
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underlying Yang Baxter Algebra is responsible for the integrability of these models. Yang 
Baxter Algebras are related to new mathematical structures often referred to as quantum 
groups which were introduced by Drinfeld (1986) and Jimbo (1985). On the other hand, 
these quantum groups have attracted attention as a powerful tool for studying properties 
of solvable systems. 

The isotropic XXX-Heisenberg model solved by Bethe (1931) corresponds to a ra- 
tional solution of the Yang Baxter equation (Baxter 1982) and is SU{2) symmetric. A 
deformation of the XXX model leads to an anisotropic XXZ-Heisenberg model related 
to trigonometric Yang Baxter solutions. A version of the model with open boundary 
conditions is quantum group invariant. The Bethe ansatz method was used to solve the 
quantum invariant spin chain for open boundary conditions by several authors (see e.g. 
Alcaraz et al 1987, Cherednik 1984, Sklyanin 1988, Mezincescu and Nepomechie 1991, 
Martin and Rittenberg 1992, Destri and de Vega 1992, Foerster and Karowski 1993). A 
quantum group invariant version with periodic boundary conditions has been constructed 
and analyzed in (Karowski and Zapletal 1993, 1994). We restrict here our interest to a 
chain with open boundary conditions. 

For generic values of q the representations of slq{2) are known to be equivalent to the 
ordinary SU{2) representations (Luztig 1989, Rosso 1988). However, this correspondence 
holds only if the deformation parameter q is not a root of unity {q"^ ^ ±l,r = integer). 

In contrast to this generic case the quantum group representation theory for = ±1 
is more complicated (Luztig 1989, Reshetikin and Smirnov 1989, Pasquier and Saleur 
1990, Reshetikin and Turaev 1991), because there exits two types of representations. The 
representations with non-zero q-dimension (see section |^ equation(|TB|)) are called "good" 
(Reshetikin and Smirnov 1989) or of type-II (Pasquier and Saleur 1990). Moreover, all 
of them have positive q-dimension, if and only if g = e*'^/^. They are irreducible and 
possess the same structure as the usual SU{2) ones. There are only a finite number 
of "good" representations, namely, those with spin j = 0, 1/2, 1, . . . ,r/2 — 1 if g'' = ±1. 
Those with vanishing q-dimension called "bad" or type-I representations. Some of them 
are irreducible and have spin j = [nr — l)/2. The others can be described as a mix- 
ing of two representations of the generic case. They are reducible but indecomposable. 
This phenomenon has a consequence for the eigenstates of an quantum group invariant 
Hamiltonian. If q approaches a root of unity, some eigenstates which correspond to "bad" 
representations become dependent and will mix. Due to the coincidence of two originally 
independent eigenstates the whole eigenspace is not complete, i.e. the Hamiltonian may 
not be completely diagonalizable. This phenomenon was investigated by Bo-Yu Hou et 
al (1991). 

The existence of "good" and "bad" representations implies the decomposition of the 
state space of an Heisenberg chain with sites into "good" and "bad" states. We 
consider the state iteratively fused by the spin 1/2- representations. The "good" ones are 
characterized by the condition that all intermediate representations have also to be "good" 
ones. This means that the "good" states are described by a restricted (j < r/2 — 1) path 
on an sZg(2)-Bratteli diagram. A projection of the full Hilbert space to the subspace of 
"good" states is analog to the restriction of the solid-on-solid model (SOS) to the so-called 
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RSOS model (Andrews et al 1984), where the local height variables of the SOS model 
are restricted to the finite set (1, 2, . . . , r — 1). This provides a connection to the minimal 
models of conformal field theories related to critical phenomena of 2D systems with second 
order phase transitions. Indeed, the analysis of finite size corrections of quantum group 
invariant XXZ spin chains (Hamer et al 1987) and RSOS-models (Karowski 1988) lead to 
conformal charges smaller than one. 

In section ^ we present the model in terms of Pauli matrices and Temperly-Lieb 
operators. In addition we write some sZg(2)-formulas which will be used later. In section 
^ we formulate the model in terms of the path basis. This formulation is closely related 
to the RSOS model. It leads to an explicit "quantum group reduction". 

In section ^ we investigate the model by means of the Bethe ansatz method. The 
eigenstates and eigenvalues of the Hamiltonian are described by sets of parameters, the 
Bethe ansatz roots. Such set of roots satisfy a system of algebraic equations, the Bethe 
ansatz equations (BAE). As a generalization of Bethe's (1931) work Takahashi (1971, 
1972) introduced for the XXZ-Heisenberg model the general string picture. This means 
that a solution of the BAE consists of a series of strings in the form 

X = A + im, A = real, m = -M, -M + 1, . . . , M, M = 0, ±1/2, ±1, ±3/2, . . . 

with a positive parity or strings with a negative parity 

A = A±zr/2±zm. A = real, m = -M, -M±l, . . . , M, M = 0, ±1/2, ±1, ±3/2, . . . 

In this context the question arises how "good" and "bad" eigenstates are determined by 
the solutions of the Bethe ansatz equations. We conjecture that "good" states are exactly 
given by strings of positive parity restricted to the first periodicity strip < r/2 and 
whose total spin fulfill j < r /2 — 1. We have no rigorous proof for this conjecture, but 
we are able to show the following coincidence. The total number of all possible string 
configurations of restricted positive parity coincides with the number of "good" states 
in the path picture of the model, which is counted by the number of all "good" path's, 
i.e. all restricted path's on the spin 1/2 s/g(2)-Bratteli diagram. This completeness will 
be shown in section ^. Such "proofs of completeness" of Bethe ansatz solutions are well 
known for models with group symmetry. Already Bethe applied this procedure to the 
XXX-Heisenberg model. For other models see EBler et al (1992a) for the SU{2) x SU{2) 
symmetric Hubbard model and Foerster and Karowski (1992, 1993) for the spl{2, l)-t-J- 
model. After having completed this paper we received a preprint (Kirillov and Liskova 
1994) in which a "completeness proof" is treated for XXZ spin chains with periodic bound- 
ary conditions which have no quantum group symmetry. In contrast to these models, we 
here present for the first time a "completeness proof" for an anisotropic quantum group 
symmetric (spin 1/2) model where the number of "good" states is smaller than 2^. 
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2 Tensor basis description of H 

The slq{2) invariant XXZ Hamiltonian with open boundary condition can be expressed 
in terms of Pauh matrices 



N-l / I -1 



g-g 



-1 



1^1 - (^n) 



i=l 



This expression may be rewritten in terms of Temperly Lieb operators 



which satisfy the relations 
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The Hamiltonian H reads 

Af-l 

h = -2y: n. (4) 

k=l 

This formula suggests the notation of open boundary conditions. 

The model is quantum group symmetric, since the Hamiltonian H commutes with the 
generators S"^, S"^ of the q-deformed algebra f/g(s/(2)) 



0. 



(5) 



We list here some formulas which will be used later. The generators possess the following 
properties 

S-] = [2S%, q'-S^q-'' = q^'S^. (6) 
The q- number [x]q is defined as 



X 



q-q 



-1 ■ 



(7) 



In the limit g — > 1 one finds [x]g = x and @ tends to the usual relations of SU{2). For g 
equal to a root of unity one has in addition 



for g^ = ±1. 



The underlying Hopf algebra structure (Drinfeld 1986, Jimbo 1985) is described by co- 
product, antipode and counit 



A(g±^^ 
e{q^'' 



±5^ 



q 



A(gS±) = qS^ + ® q-'^' 

l{q'") = -q^'S^, 
e{q^^) = 0. 



(9) 
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The configuration space (C^)'^^ of a quantum spin 1/2 chain with N lattice points has 
the natural tensor basis 

\a) = \aN: . . . , ai) = {on) ■ ■ ■ ® \ai) 

with I — |) = (i) = (i) (spin down) and ||) = (t) = (^q) (spin up). The representation of 
the quantum group generators on the configuration space is given by 







■®q^'/^ 




(10) 
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The q-deformed Casimir operator is 

^2 = S-S+ + [5^ + 1/2]J - [1/2];. (13) 

For generic values of q all representations pj arc irreducible and highest weight ones with 
spin j = 0, 1/2, 1, . . . (see e.g. Lusztig 1989, Rosso 1988). If q tends to a root of unity 
(g** = ±1) the Casimir eigenvalue 

s'-\j + ml - ml (14) 

can take identical values for different spins j and / in case 

j' = j + nr or j' = r — 1— j + nr, n E Z. (15) 

The representations pj and pjr, which are different for generic q, mix for q^ = ±1 and 
build reducible but indecomposable representations. Together with the irreducible rep- 
resentations with spin j = {nr — l)/2 they are characterized by vanishing q-dimension 
defined by 

dim,p = try,(g-2^') (16) 

where is the representation space. Therefore they are called "bad" or type I repre- 
sentation. On the other hand, there exits a finite number of representations which are 
in one-to-one correspondence to SU{2) representations. They are irreducible, have spin 
j = 0,l/2,...,r/2 — 1 and non- vanishing q-dimension 

dj = dim,p,- = [2j + l]g. (17) 

They are called "good" or type-II representations. Their q-dimension is always positive 
if g = e*^/''. 

Reshetikhin and Turaev (1991) have shown that for the "good" and the "bad" repre- 
sentation spaces the following structure holds 

Vgood^Vgood = (0V^pood)e (0Had), (18) 
Vgood^Vbad = 0Had- (19) 
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This means that "unitarity" in the "good" subspace (see also Reshetikin and Smirnov 
1989) is fulfilled for "good" covariant operators 

( good' I I good ) = ^ ( good' \ Aj^^^^ \ good" ){ good" \ Bk^^^^ \ good ). 

good " 

(20) 

In the next section the "good" representations in the state space (C^)®^ of the XXZ- 
Heisenberg model are characterized in terms of the "path picture" . For q equal to a root 
of unity the state space may be reduced to the subspace of all "good" states. Changing 
the metric in this subspace the Hamiltonian (|l]) will become selfadjoint. 



3 Path basis formulation 



For generic values of the deformation parameter q the representations of slg{2) are irre- 
ducible and classified as those of the undeformed group SU{2). The space of the spin 
j representation is spanned by a set of basis vectors | j, m ), m = j, j — 1, . . . , —j with 



S'^\j,m) = ^J[J Tm]g[j ±m + l]g\ j,m±l ) 
\ j,m ) = m \ j,m ). 



(21) 



The tensor product space V^^ ® V^^ decomposes into a direct sum of irreducible spaces 





i=\h-h\ 



given by Clebsch-Gordan coefficients (see e.g. Kirillov and Reshetikin 1989) 



mi ,m2 



J 32 Jl 

m 1712 ^1 



(22) 



By successive fusion of the spin 1/2 states associated to the lattice sites we construct 
the state 

\ 3_^T^) = i,m) (23) 

a 

which is labeled by the "path" of spins j = {Jn, Jn-i, ■ ■ ■ , j2, ji = 1/2) and the mag- 
netic quantum number m = J2iC(i- The matrix element is a product of Clebsch-Gordan 
coefficients 



( a I j, m ) = 

m2,...,mpf_i 



Jn 1/2 Jn-i 



J3 1/2 J2 

ms as 1712 



32 1/2 1/2 

m2 012 0L\ 



The spins jfc are restricted by the fusion rule 

jfc+i = jfc ± 1/2 



(24) 



(25) 
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and j = Jn describes the total spin of the state. The space is decomposed into a 
direct product of a space described by paths Wj and a space Vj where the generators 
S^,S' act 

= E ® V,. (26) 

j 

Quantum group invariant operators as the Hamiltonian only act in the path space Wj. 
By the Wigner-Eckart theorem the magnetic quantum number m is not be changed and 
in addition the path space (or reduced) matrix elements do not depend on m. Therefore 
we will omit it in the following. 

For example the Temperly Lieb operators (||) in path space act as 



Tk\ ...,jk,---) = Sj,_,j,^, J2 \ -■ (27) 

Thus, we easily obtain the matrix in path space of the Hamiltonian H, which decomposes 
into different blocks for different total spin j. 

Now we turn to the case where the deformation parameter g is a root of unity 
(g** = ±1). In the path picture it is very simple to characterize the "good" states. From 
section ^ it is obvious that the "good" states are given by all restricted paths 

U3ood)=l-^----'-^i)' 2jfc + l<r, k = l,...,N. (28) 

For generic values of q the number of all states with total spin j is equal to the number 
of all possible unrestricted paths (p5|) on the s/(2)-Bratteli diagram 

^^={n/2-j)~{n/2 + 1+j)- ^^^^ 

For = ±1 the number of all "good" states with total spin j < r/2 — 1 is equal to the 
number of all possible restricted paths ( ^81) on the cut s/g(2)-Bratteli diagram 

oo 

Qj = E ^j+rk- (30) 

fc=— oo 

In the next section, it turns out that this number coincides with the number of, what we 
will introduce, the "good" Bethe ansatz states. Note that only a finite number of terms 
contribute, because the binomial coefficient is defined to be zero for integer n < or 
n > m. 

The Hamiltonian does not lead to a transition from a "good" state to a "bad" one. 



This follows from relation (pTf). The only possible "good" path's | j ^^^^ ) which lead by 
an action of H = —2 J^Tk to a. "bad" state | j'^^^ ) by 

I l,oo, ) = ^ I l,oo, ) + ^' I ) (31) 
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must have one or more values jk+i = jk-i = jfc + 1/2 with 2jfc+i + 1 = r — 1. Then 
the "bad" path would have the value of = Jfc + 1 = ('^ — But the coefficient c' 

vanishes because of eq. (p7|) and dj'^ = [r]q = 0. That means, the matrix H in path space 
decomposes into a "good" submatrix and a "bad" part. Note that these arguments cannot 
be directly translated to the tensor picture, because the transition matrix elements (p^ ) 
become singular for "bad" states. 

We now introduce a new scalar product in the "good" subspace in the path formulation, 
such that the Hamiltonian becomes self adjoint 

{ l,rn' \ i,m) = Sm'mSfj_. (32) 

Note that this scalar product coincides with the natural one in the tensor picture only for 
real values of q. The selfadjointness of the Hamiltonian follows from eq. (|27| ) 

= H for q = e*"/^ r = 2, 3, . . . (33) 

since the q-dimensions dj = [2j + l]g are positive for 2 j + 1 < r. 



4 Bethe ansatz method 

The eigenstates and eigenvalues of the XXZ-Hamiltonian ([l|) with open boundary condi- 
tions are described by sets of distinct spectral parameters {Ai, . . . , A;}, the roots of the 
Bethe ansatz equations (BAE) 

/^sinh7(Aj +i/2)y^ _ ' sinh7(Aj - Afc + i) sinh7(Aj + Afc + i) 



\sinh 7(Aj - i/2) J fc=i,Mi ^^^^ ^(^i ~ - i) sinh -f{Xj + A^ - i) 

where q = exp(i7). Because of the open boundary conditions it is sufficient to consider 
only roots with positive real parts < JJAfc < oo. The energy is 

sin^ 7 , , 

E = -4y . 35 

cosh27Afe - cos 7 

The total spin j of an eigenstate is related to the number of roots / by 

j = N/2-1, 1 = 0,..., N/2 (36) 

where the lattice length is assumed to be even. Because of the quantum group invari- 
ance the Bethe ansatz solutions are highest weight states (Destri and de Vega 1992). The 
aim of our investigation is to find a concrete criterion for the Bethe ansatz solutions to 
be "good" states in the sense of section 2 and 3. 

A detailed numerical analysis of the BAE motivated also by the observation of Destri 
and de Vega (1992) that when q approaches a root of unity the appearance of a "bad" 
state correspond to a root tending to infinity (see also Appendices A and B) leads us to 
the following 

Conjecture 1 : For q = exp(z7r/r) a Bethe ansatz state is the highest weight vector of a 
"good" representation (in the sense of section^ and^ see eq. /jl^j, if and only if 
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(i) the total spin j is restricted by 2j + 1 < r, i.e. the number of roots I must be larger 
than {N +l~r)/2, 

(a) the roots are restricted to the first periodicity strip \ 53 |< r/2. 

The first condition (i) is obvious from tlie definition of "good" states. We liave no 
rigorous proof for the second condition (ii), but we can show completeness in the sense 
that the number of "good" Bethe ansatz solutions defined by Conjecture 1 leads to the 
correct number of all "good" states on the lattice counted by all paths on the Bratteli 
diagram (see eq. (|30D). To "prove" Conjecture 1 we proceed as follows. We classify 
the Bethe ansatz roots by means of Takahashi's string picture. A configuration of strings 
is given by sets of integers. In Conjecture 2 we give upper bounds for these integers. 
This is also for the group symmetric case a nontrivial (but simpler) problem. The upper 
bounds are smaller than a naive estimate would suspect. The reason for this phenomenon 
is that the string picture is only a very rough approximation. The exact set of roots are 
given by deformed (sometimes even degenerated) strings. This leads to the fact that less 
roots are possible than the exact string picture would allow. For the group symmetric case 
already Bethe (1931) was able solve this problem, because there is a natural way to fix 
the bounds in order to get the correct number of states. This is much more complicated 
for the quantum group case. Only guided by a lot of numerical calculations, we were able 
to solve this problem. 

It is convenient to classify solutions by the so-called string hypothesis (Takahashi 1971, 
1972). Any solution of the BAE consists approximatively of a series of strings in the form 

AM = AM + im, m = -M,-M + l,...,M, M = 0, ±1/2, ±1, ±3/2, . . . (37) 

with positive parity or strings with negative parity A = A±ir/2± im. According to 
Conjecture 1 (ii) we must take into account only the positive type (|37|). 
A Bethe vector is characterized by the sets of real string centers 

{{Ao}, {Ai}, . . .} where {Am} = {A^.i, • • • , ^m,um}- (38) 

The length of a string is 2M + 1 and i>m is the number of (2M ± l)-strings with different 
centers A^./c- The total number of roots writes as 

l=Y.{2M + l)z/M. (39) 

M=0 

As usual we rewrite the Bethe ansatz equations in terms of the strings centers 

2M 

n k^(2Am.) 

m=l 

^ n n yM,M'i^M,t - ^M',k) VM,M'i^M,i + ^M',k), (40) 
M' ^=1 

{M,i}7^{A/',fc} 
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where 

Vm,m'W= n Kn(A)Kn+i(A) and K.(A) = . , ^ • (41) 

Taking the logarithm of eq. ( ^0]) sets of integers {Qm} occur which determine the solutions 

2M 

2A^^M+l/2(AA/,^) = 27rQM,, + 2^'„(2Am,^) 

m=l 

+ E E (*Af,M'(AA/,*-AA/',fc) + ^M,M'(AAf,i+AAf',fc)j (42) 



M' ^=1 

{A/,i}#{A/',fc} 



where 



Af+A/' 

^Aw(A) = (^m(A) + ^,„+i(A)) and (43) 

m=|Af-A/' 

^m(A) = 2arctan(cot(7m)tanh(7A)) , m > 0, ^o(A) = 0. (44) 
A solution of the BAE is parameterized by a configuration of distinct integers 

< Qa/,1 < Qm,2 <■■■< Qm,um < Q'^r^ ^ = 0, 1/2, 1, . . . (45) 



max 
A/ 



which can be used to count the total number of possible solutions. The upper bounds Q 
are defined by the following condition. All roots corresponding to the integers fulfilling 
(1^) are "good" ones in the sense of Conjecture 1. However, if one Qm,um — Ql^F^ + 1 
then one root would be pushed to oo or it would return as a "bad" one with \'^\\ > r/2. 
This is meant in the sense of deforming the model with respect to q. The counting of the 
"good" states may be performed by the following 

Coniecture 2 : For "good" Bethe ansatz solutions: 
(in) the upper hounds in relation are given by 

Qtr = '2j + i^M+2 E 2iM'-M)uM'-GM, Ga/ = max(2j+2M+3-r, 0) (46) 

M'>M 

(iv) and the string length is restricted by 

2M + 1 < r - 2. (47) 



Note that this condition (iv) is stronger than expected from Conjecture 1 (ii) and 
which would allow also strings of length r — 1. However, (iv) follows from (iii), since for 
the maximal string the sum in eq. (H6[) is empty and by relation (|45D 
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Note also that a naive estimate would lead to a larger upper bound than that of (iii), 
namely to the largest integer smaller than Q^, which is obtained by putting A^f^j oo 
in eq. (^) (using \l/m(oo) = vr — 2jm) 

= 2j + z/M + 2 ^ 2(M' - M)z/M' + (2M + 1) (l - -{N - 2/ + 2M + 2)] . (48) 
M'>M \ r J 

The upper bound of (iii) is obviously smaller than that obtained from this equation, even 
for the isotropic case g = 1 or r = oo where Q™l^{q = 1) = Qfj — 2M — 1, already used 
by Bethe (1931) with J{M,M') = 2min(M,M') + 1 - 5m,m'1/2 and 

QM'(g = l) = 2J + Z/M + 2 J2 2{M'-M)iyM' (49) 

M'>M 

= N -2J2''M'J{M,M'). (50) 

M' 

We are not able to prove (iii) rigorously, however, we performed many numerical calcula- 
tions (see Appendices A and B) and justify it by the following counting of "good" states 
leading to the correct result. 

We must solve the following combinatorial problem to compute the number of Bethe 
ansatz states. The integer Q'^^ denotes the number of vacancies for a string of given 
length 2M + 1 . The number of possible configurations ( ^5]) is {^i^^^ ) and therefore the 
number of combinations of a given set {um} reads 



imax \ 



z(iv,{M,r)=nr'' )• (51) 



M 



Now we obtain the total number of Bethe ansatz states with fixed / by taking the sum 
over all configurations {pm} 

Z{N,l,r) = Z{N,{uM}.r) with fixed / = ^(2M + 1)z/m. (52) 

{fM} M 

Extending the calculations made by Bethe (1931) one can show for M > 

gM^(iV,{^M},r)=Q7^i/2(A^-2/x,{//M},r-l), /x = E ^Af, i^'m = (53) 

M 

This implies with eqs. ( ^T]) and (|46|) 

Z(Ar,{z/,,},r)= n° )z(Ar-2;U,{z/;,},r-l) with Q^^^ = N -2^i+u^-G,{r). {U) 

We introduce the partial number of configurations depending on the number of strings /i 

Z{N,l,ii,r)= Z{N,{pM}.r). (55) 

^(2M+l)vM='- 
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From eq. we obtain the following recurrence relation 



Z(iV,/,/i,r)= [ )^(iV-2/x,/-^,/x-z/o,r-l). (56) 



uo=0 



1^0 



This relation differs from that of Bethe for the XXX-model by the additional dependence 
on r. Note that the equation (EB) holds for fi < I. 



The initial values of this recurrence relations are given hj I = fi where only real roots 
exist {uq = I and i^m = for M > 0), therefore from eqs. ( ^T] ) and 



ZiN,l,,^l,r) ^ + (57) 

Z{N,l,fx,r = 2) = 0. (58) 



The second initial condition for r = 2 follows from eq. 

We now introduce the functions fk,d depending on the integers k and d = 0,1 

( N-l-kir-2)-Goir + l) + d \ /l + k{r - 2) ~ d\ 
Lemma : The function 

CO 

Z{N,l,f^,r)= Y (/m-/m) (60) 



solves the recurrence relation (|3^ and fulfills the initial condition ([3^. Thus it is equal 
to the partial number of states defined in eq. ([5^j. 



The proof of this lemma is performed in Appendix C. 

In the isotropic limit r ^ oo all the functions fk^i — fk,Q vanish except that for A; = 
which coincides with the solution already found by Bethe (1931). One can calculate the 
total number of Bethe ansatz states for a given number of roots / (E^) by 



Z{N,l,r) = j2Z{N,l,fi,r). (61) 
For the interesting case 2j + 1 < r (with j = N/2 — I) it has the form 

ZiN,l,r)= Y r,+,., r,+,., = Y: {fk,i - fk,o) = , _ . J - / _ 1 - J (62) 

k=—oo fj.=l \ / \ / 

which indeed coincides with the number of "good" path's (0). Thus, we have shown that 
the conjecture about Bethe ansatz states of the "good" type yields the correct number of 
states with respect to the configuration space of the quantum invariant spin chain restricted 
to the subspace of "good" representations . 
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5 Conclusions 



The configuration space of tlie quantum group invariant XXZ Heisenberg model witli open 
boundary conditions decomposes into a "good" and a "bad" part, if and only if the de- 
formation parameter g is a root of unity. Furthermore, if q takes the values q = exp(i7r/r) 
(r = 3,4,...) a new metric may be introduced in the "good" subspace such that the 
Hamiltonian becomes selfadjoint. This is reminiscent of minimal models of conformal 
field theory. The central charge of the Virasoro algebra is known to be restricted by 
unitarity (Friedan et al 1984) to the values 

c= l-6/r(r-l), r = 3,4,5,.... (63) 

Indeed, by finite size computations (Hamer et al 1987) of the spin chain one obtains the 
identification to c. 

In this paper the completeness of the Bethe ansatz of the quantum group symmetric 
spin 1/2 Heisenberg model was proved for a configuration space reduced to the "good" 
representations at g = exp(i7r/r) with r = 3,4,5,.... The "good" Bethe ansatz solu- 



tions (in terms of strings of positive parity) are parameterized by integers (|45| ) which 
are bounded by the upper values Q™/'' (^61) . The conjecture for Q^^^ in the quantum 
group case introduced in this paper yields the correct number of states. Furthermore, the 
"good" Bethe ansatz solutions were checked numerically for small lattice length (A^ < 12) 
by solving the BAE exactly. In addition the eigenvalues obtained by these solutions were 
compared with those obtained by diagonalization of the Hamiltonian in the path basis 
(see eq. (P). 



Appendix A: Non-string solutions of the BAE 

The Bethe ansatz equations in the logarithmic form (^) provide a determination of 
solutions by a set of integers {Qm} if one assumes that complex roots consist of strings 
A = A±im with known imaginary part m. However, this assumption is known to be only 
an approximation. In general we have complex pairs X = x±iy where the imaginary part 
is not a integer or a half-integer (in case of finite lattice length A^). To consider exact 
solutions we analyze the BAE (0) without this conjecture. Unfortunately, for / > 2 (/ . . . 
number of roots) such an investigation can be made only numerically. But analyzing a 
large amount of numerical solutions we were lead to some rules which should be viewed 
as general properties of the Bethe ansatz equations. 

We have found (see appendix B) that by a successive increasing of 7 to integer values 
r = 7r/7 (g = exp(i7)) the largest admissible value of Q must be reduced in order to 
obtain finite rapidities of positive parity (which are assumed to be related to "good" 
states). Especially, if a 2-string {I = 2) tends to infinity at Q and r then a 2-string at 
Q' = Q — 1 degenerate at r' = r — 1 into two 1-strings where one 1-string tends to infinity. 
This leads to the corrections G1/2 and Go of the upper boundary Q™^^. In the general 
case / > 2, we can also assume that the largest value Q^^^ for a set of integers is caused by 
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non-string effects. By numerical computations we have observed the following important 
property. 



If a k-string (string of length k=2M-|-l) tends to infinity at Qm and r then a 
k-string at Q'j^j = Qm — 1 and r' = r — 1 degenerates into a divergent (k-1)- 
string and a 1-string. Moreover, a k-string at = Qa/ ~ 2 and r" = r — 2 
degenerates into a divergent (k-2)-string and a 2-string. In general, a k-string 
at Q'm = Qm — n and r' = r — n degenerates and the new (k-n)-string or the 
n-string becomes infinite. 

Due to this conjecture we obtain the condition 

GM{r)=GM{r-l)-l. (64) 
Considering the value Q^ in equation (|^) 

Q^j = N-2y2 ^M'J{M, M') + (2M + 1)(1 - ^(A^ - 2/ + 2M + 2)) 
M' ^ 

one can assume that the initial value Gm{t) = 1 is given by 

-(iV-2/ + 2M + 2) = 1 (65) 

TT 

which leads to 

GM(r) =2j + (2M + l)-r + 2 (66) 



as was suggested in ([461). 

We point out that this picture of admissible strings holds only in the case 2j ' + 1 < r. 
For 2j ' + 1 > r there are more configurations which leads to finite solutions of positive 
parity. Nevertheless, such states always belong to the "bad" type and have to be excluded. 

Appendix B: Numerical computation of the BAE 

Considering the two-magnon case / = 2 one obtains two kinds of possible solutions. Due 
to (|3^ ) there are either two distinct real roots {Ai, A2} or a complex pair {x + iy,x — iy}. 



The two real roots satisfy the following Bethe ansatz equations 

2'kQi = 2Ar^i/2(Ai)-^i(Ai-A2)-^i(Ai + A2) (67) 
27rg2 = 2A^^i/2(A2) -^i(A2- Ai) -^i(A2 + Ai) (68) 

corresponding to two 1-strings. Because of the zero imaginary part the solution is exactly 
given by (^2|) with z/q = 2 and i/^ = (A; > 0). If one root tends to infinity the value Q 
reads 

goo^2j + 2-^(2j + 2). (69) 

TT 



14 



Thus, for {2j + 2)71/7 < 1 the largest admissible integer is given by 

g-- = 2j + 2 (70) 

and for Tc/'-f = r = (2j + 2) one root tends to infinity. The condition 

Q-- = 2j + 2 - Go, Go = 1 = 2j + 3 - r (71) 

ensures finite solutions. For example, the case = 8 (/ = 2, j = 2) and r > 2j + 2 
provides configurations of integers {Qi,Q2} which are in one-to-one correspondence to 
the solutions of the BAE 

{Ai,A2}^{Qi,Q2} 

with 1 < Qi < Q2 ^ Q'q^^ = 2j + 2 = 6. At r = 2j + 2 the number of configurations is 
reduced by Q™'^^ = 6 — 1. For r < 2j + 2 we have no "good" state. 

The case of complex roots (Z = 2) is more complicated. Using the string conjecture 
we would have the following equation for the complex pair 



7rQi/2 = Ar^i(Ai/2) - ^i(2Ai/2) (72) 

with the configuration z/1/2 = 1, z/o = 0, z/i = . . .. As a function of A the r.h.s. is 
monotone and all integers Q up to the asymptotic value 

Q^2 = 2j + z/i/2 + 2-2(2j + 3)/r 

would lead to a finite value Ai/2- This contradicts the condition that the largest admissible 
integer is given by 

Q^/T = 2j + iyi/2 - Gi/2 = 2j + 1 - Gi/2 (73) 

Even in the XXX case (G1/2 = 0) there are integers Q with Q^yT < Q < QT/2 which 
are allowed in (|T^ but forbidden as Bethe ansatz states. We have already pointed out 



that for such large values Q the imaginary part of the root tends away from the assumed 
(half) integer value and therefore the string conjecture does not hold shown by the exact 
solutions of the BAE. A complex root denoted by A = x + iy satisfies the equations ( ^4]) 



/ sinh7(x + i(±?/ + l/2))\^^ sin7(±2y + 1) sinh7(2a; -|- i) 
I sinh7(x + i(±2/ — 1/2)) j sin7(±2y — 1) sinh7(2a; — i) 



(74) 



Considering the phase of these relations 



vrQi/2 = -N<Pi/2±y{x)- N(Pi/2^y{x)+<Pi{2x) + Tx{N -I) (75) 
'Pyi.^) = 2 arctan(tan 7?/ coth 7x) (76) 

we obtain in the limit y —>■ 1/2 (with (pyix) = tt — ^j^(a;)) the equation for a two-string 
(ff^). Thus, these integers Q are related to (|T2|). The magnitude reads 

(77) 



/ sinh^ 7X + sinh^ 7(±1/ + 1/2)) \ sin^ 7(±2y + 1) 



\sinh 7X + sinh ■j{±y — 1/2)) J sin 7(±2?/ 
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which is used to ehminate the real part x in (|75[). With 

X = b{y) {7i 



wUy) = + (80) 

sinh 'ymyy — 1/2) 
the counting function z{ii) depending on y reads 

7,z{y) = -N (0i/2+,(%)) + 0i/2-,(&(2/))) + 0i(2K?/)) + <N - 1). (81) 

It turns out that a solution for | y |< 1/2 and | y |> 1/2 is given by even and odd integers, 
respectively 

Qi/2 = z{y). (82) 




Figure 1: Counting function for a complex pair A = x ± iy with N = 8, I = 2 and r = 14 
as a function of y. Full (open) circles correspond to odd (even) values of Q. 



The general behavior of z drawn in figure |^ shows that for small values of Q the 
imaginary part y of a complex pair is approximated very well hj y = 1/2 with the 
correspondence 

{x + iy,x- iy} ^ {Qi/2}- 

But for large values Q the imaginary part of the root tends away from the line y = 1/2. 
Now we discuss the example N = 8 but it can be easily generalized. The counting 
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function z{y) takes it maximum at y = (for y < 1/2) and at y = ir/A'y = r/A. One can 
show that for 7—^0 (XXX case) the maximal integer is determined by (5™2^ = 2j + 1 
which corresponds to the upper restriction stated above. Ahhough the maximum z™^^ 
(at y = r/4) is larger than 6 

6 < z'"'^" < 7, y> 1/2 
the largest admissible integer is related to 

= 2j + 1 = 5. 

Note that the odd integers Q1/2 = 1, 3, 5 determine the solutions with y > 1/2. The even 
integers Q1/2 = 2,4 correspond to roots with y < 1/2 because of 

4 < < 5, y< 1/2. 

The maximum of z{y) for both y < 1/2 and y > 1/2 depends on the anisotropy 7. It 
decreases with increasing 7. The picture described above remains unchanged up to values 
r > 2j + 3. In the case r = 8 for example we have z™^^2 = 5.25 and ^^1/2 ~ 4.07195 . . .. 
Thus, the integers Q1/2 = 1, 2, 3, 4, 5 are assumed to correspond to "good" state if r > 8. 
Increasing 7 we have at 

7r/7 = r = 2j + 3 = 7 

(according to (pT))) a divergent complex pair with Q1/2 = 2j + 1 = 5 {x 00, y = ±r/4) 
because of -2^1/2 ~ ^- Therefore, the value Q™T finite roots must be reduced by one 

= 2j + 1 - 1 = 4 

which leads to 6*1/2 = 1- The odd value Q1/2 = 2j + 1 = 5 is now forbidden. Furthermore, 
by increasing 7 or increasing N we reach a critical value where the maximum of z{y) for 
y < 1/2 is smaller than the even integer Q1/2 = 2j = 4 

<<V2 = 3.93162... <2j, r = 7. 

The imaginary part y tends to zero. This fact was observed by Vladimirov (1984) and 
EBler et al (1991b) in the XXX case. The complex pair was found to be replaced by two 
additional real roots. The dependence on 7 was investigated by Jiittner et al (1993). Here 
we have a similar situation. The additional real solution is considered as a degenerate 
complex pair corresponding to Q1/2 = 2j = 4. The roots are finite and therefore remain 
"good" . Such a degenerate solution is described by (|68D having identical integers 

Qi = Q2 = Qr' = 2j + 2 = 6 

at Q^^^ although their rapidities are distinct Ai 7^ A2. This correspondence is denoted 
by 

{Al,A2}--{Q^^Q^n--{Ql/2}. 
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Note that there is a difference between Q™^^ of real roots and the value Q1/2 for a 2-string. 
At r = 7 we consider the values Q1/2 = 1, 2, 3, 4 as "good" integers. 

Increasing 7 to r = 6 we will see that a solution related to Q1/2 = 4 becomes "bad". 
It can be shown that in the limit 71/7 — >■ 2j + 2 a complex pair at Q1/2 = '^j = 4: always 
degenerate. But at r = 2j + 2 = 6 a real root related to Q = 2j + 2 = Q™^^ = 6 is known 
to tend to infinity. This happens also for a degenerate complex pair. Because one part 
of this degenerate string is divergent we must exclude the corresponding integer of the 
2-string 

Q^;^ = 2j-i 

which leads to G1/2 = 2 and "good" values Q1/2 = 1, 2, 3. 

A further increasing of 7 to r < 5 leads to the case where any state belongs to the 
"bad" sector. Namely, the total spin (j = 2) is related to 2j + 1 > r. One can check that 
at r = 5 the admissible numbers for finite solutions of positive parity read Q1/2 = 1, 2, 3. 
Thus we have no further reduction of Q™^^- -^^^ ^^^^ '^^ interest for our consideration 
of "good" states. 

The case / > 2 is treated as follows. Fixing a given set of integers {Qm} we first 
calculate the centers of strings {AM}-yo for < 70 << 1 by fixed point iterations of the 
BAE (H) 

if / 2M i^M' \ \ 

AM,i = - tanh"^ tan 7m tan 27rQM,i + E 2^™(2AM,i) + E E • • • /^^ (83) 

7 \ \ m=l M'=Ok=l ) ) 

It turns out that the iteration converges (nearly) independently on the initial values for 
AM,i- The set of centers {Am}7o is now a function of {Qm}- The next step consists in 
an exact computation of the roots {A} by the original BAE (0) which is numerically 
solved by the Newton method. As initial values we use the string centers {Ajv/j^g where 
the corresponding imaginary part of each member reads 

A = Km + im + i(5 

with <| (5 |<< 1. With a careful choice of 5 the Newton method converges which yields 
the accurate solutions of the BAE for 79. Thus, we have a correspondence between the 
set {Qa/} and {A} 

{^}7o = /({<5a/},7o)- 
Now the anisotropy 7 is increased in small steps 

7fc+i = 7fc + A. 

For each 7^ the Newton method is applied. In contrast to the first step we now take initial 
values which are accurate solutions in the step k — \. If A is small enough one can assume 
that the roots {Aj^^. keep the relation to {A}^j._^ with respect to the correspondence to 
the integers {Qm}- This provides 

{A},, = /({Qm},7.) 
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Table 1: Exact solutions for N = 8, 1 = 3. 





= 1 


= 2 


= 3 


71-/7 


Ai 


^2,3 


Ai 


A2,3 


Ai 


A2,3 


7 


0.672828 


0.672829 ±il. 000377 


1.305862 


1.338141 ±il. 035500 


2.379542 


2.183152 ±il. 348160 


6 


0.716146 


0.716398 ±il. 000649 


1.407535 


1.489400 ±il. 060994 


oo 


oo ± iir/Sj 


5 


0.810005 


0.811995 ±il. 001089 


1.681845 


oo ± z7r/47 






4 


oo 


1.376812 ±i0.569724 











which can be used to investigate the behavior of the solutions as a function of the integers 
{Qm}. 

Now we discuss some numerical examples. Considering a 3-string for = 8 and / = 3 
{ui = 1) the integers Q = 1,2,3 are allowed {n/'y > 6) and the exact solutions are listed 
in table |l|. The numbers Q = 1,2,3 are associated to the following picture where the real 
and imaginary part of the roots are plotted in the x-y-plane. At 7 = vr/? we obtain 



1 




_ * 



Q=l 



Q=2 



* 

Q=3 



The roots are arranged as 3-strings very well. Increasing 7 to tt/Q the roots related to 
Q = 3 tend to infinity - the real root as well as the real part of the complex pair. 



1 




Q=l 



Q=2 



Q=3 



Thus, the admissible numbers are reduced to Q = 1,2. Now, in the limit 7 vr/S the 
string at Q = 2 degenerate. The real root remains finite - but the difference to the real 
part of the complex pair becomes larger. This 3-string degenerate into a 2-string and a 
1-string in such a manner that at vr/S the complex pair tends to infinity. 



1 




Q=l 



Q=2 



"bad" 



Q=3 
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Now we have only one allowed integer Q — 1. 

The behavior of a 4-string is similar. If N 
admissible values Q — 1,2,3 {j < it/7). 



— 10 and I — 4, i/3/2 = 1 we have 3 



1 




* 
* 



* 
* 



* 
* 



Q=l Q=2 Q=3 

At 7 = 7r/7 the configuration related to Q = 3 becomes infinite 



1 





* 



Q=l 



Q=2 



Q=3 



and the allowed integers are reduced to Q = 1,2. If the anisotropy is increased to a 
critical value between 7r/7 > 7 > 7r/6, the imaginary part of one complex pair becomes 
smaller and vanishes at the critical value. This complex pair is replaced by a pair of two 
real roots. 

* 

1 * 







X 



Q=l Q=2 Q=3 

Thus, the 4-string degenerates. The whole configuration can be viewed as a set of one 
finite 1-string and one 3-string which tends to infinity at 7 ^ 7r/6. 



1 




"bad" 



Q=l 



Q=2 



Q=3 
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In this case only Q = 1 is permitted. However, in the hmit 7 — 7r/5 this set degenerates 
into two (2-string hke) complex pairs. One 2-string is divergent. 



1 




* 



"bad" 



'bad" 



Q=l 



Q=2 



Q=3 



Therefore, no integer is allowed for a 4-string (71/7 = r < 5). 

These examples and similar ones lead us to the conjectures about admissible configu- 
rations of finite strings of positive parity stated in section H. 



Appendix C: Proof of the Lemma in Section g| 

We define the function 

00 

Z(iV,Z,/i,r)= J2 (/m-/m) 
fc=— 00 



and show that it solves the recurrence relation (|56|) and fulfills the initial condition (|57D . 
Only a finite number of terms contribute because the binomial coefficients (^^^ are zero 
for m < or n < 0. 

The right hand side of the recurrence relation reads for ^ 



2^ ( l/fc,d(A^-2/i,/-/i,/i-z/o,r-l) (85) 



1^0=0 



g fN-2^ + uo-Goir) 



uo=0 



' N -I - k{r -3) -Go{r) + d- fi ' 
.N-l + l-k{r-2)-2fi- Go{r) + z/q , 



n + k{r -3) -d- fi^ 
+ k{r - 2) - 2/i + z/oy 

Using the sum rule of Binomial coefficients 



^6) 



'a + b^ 
. n J 



E 

fe=0 



,k \n — k I 



^7) 



one has (with g = Go{r) and u = z/q) 



E 



'N -2i2 + u -I - g + 1 - k{r -2)\fl- 1 + k{r -2y 
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which is inserted in Fk^d- With the identity 

^N-2f^ + u-l-g + l-k{r-2)\^ N-l-g + d-k{r-3)-^i ^ 
\ u-uj ) \N - I + 1 ~ g - k{r - 2) - 2fx + 

_ /N-l-g + d-k{r-3)-fi\fd+k + fi-l-u\ 

\ d + k + ^— 1 — LU J \ U — LU J 



^9) 



we have 



Fk,d — 



oj=0 



'l-l + k{r-2)\/N-l + d-g-k{r-3)-iJ> 
. u I \ d+k + fi — 1— uj J 



^ fd + k + fi - 1 - uj\ f I - d + k{r - 3) - fi\ 

k.[ )[l + kir-2)-2f, + u)- ^ ' 

With the help of the sum rule (|87D this equation reads 

'I - 1 + k{r - 2)\ /N - I + d - g - k{r - 3) - fi^ 
^ LU J \ d + k + fi — 1 — uj y 

' l-l + k{r-2)-tu \ 
- 1 + k{r - 1) + d - jj,) 



uj=0 



X 



(91) 



^n-l + k{r-2)\/fi-d-k\/N-l + d-g-k{r-3)-fi\ ^g^) 



1^=0 \ 

with the final result 



H — d — k I \ CO I \ d + k + jj— 1 — to 



' N-l-k{r-2)-g \f I - 1 + k{r - 2) \ 
- I + 1 - k{r - 1) ~ g - 12 - d) \l + k{r - 1) - fi - 1 + d) ' 



(93) 



In the "good" sector 2j + 1 < r — 1, where by ( ^U] ) Go{r) = Go{r + 1) = (7 = it is easy 
to show the recurrence relation (^) for term in (^0]) separately 

Fk,i — Fkfi = fk,i — fk,Q- (94) 

Thus the function Z{N, I, /x, r) ( P^ fulfills the recurrence relation. But in order to use the 
initial conditions (^) and ([58|) we also have to consider the sector 2j + 1 > r — 1 where 
Go{r) > and Go{r) = Go{r + 1) + 1 which yields additional terms 

Fk,i - Fkfi = fk,i - fk,o + Kk,i - Kkfi (95) 

with 

^/ N -l-k{r-2)-Go{r + l)-d \n + k{r-2) + d-l\ 

liV-/ + l-A;(r-l)-Go(r + l)- Jl I + k{r - I) - ' ^ ' 



The function Z (|8^) fulfills the recurrence relation also for 2j + 1 > r — 1 due to the 
following identity 

iKk,i - Kk,o) = 0. (97) 

k 
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We have no analytic proof of this identity but we have verified it numerically for all possible 
cases for < 60. It is easy to show that Z assume the initial conditions. Therefore, the 
function Z{N,l, fi,r) = Z{N,l, fi,r) is equal to the partial number of states defined in 
eq. (HD. 
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